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Abstract 

The problem of existence and constructing of integrals of motion in stationary quantum mechanics 
and its connection with quantum chaoticity is discussed. It is shown that the earlier suggested 
quantum chaoticity criterion characterizes destruction of initial symmetry of regular system and 
of basis quantum numbers under influence of perturbation. The convergent procedure allowing to 
construct approximate integrals of motion in the form of non-trivial combinations depending on 
operators (q,p) is suggested. Properties of the obtained integrals with complicated structure and 
the consequences of their existence for system's dynamics are discussed. The method is used for 
explicit construction and investigation of the approximate integrals in Henon-Heiles problem. 

o 

1 Introduction 

03 
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■ For last dècades the investigat ions of quantum chaos have been extensively carried out but 

this field remains under hot discussions. One part of researchers believe that quantum chaos 
^ . doesn't exist and at best should be studied in the semiclassical approximation. It is obvious 

from correspondence principle that quantum counterpart of classical system should have 
properties reflecting regularity or chaoticity of classical trajectòries. The law of level spacing 
distribution is considered to be one of such properties. It is believed that the quantum 
analogue of the chaotic classical obeys Wigner level spacing law, while Poissonian law holds 
for regular systems. However many authors (including us [1]), pointed out the incompleteness 
and crudeness of this criterion of quantum chaoticity. 

In this paper we continue to develop our approach [2-4] to the chaotic properties of the 
quantum Hamiltonian systems. Our main point is connection between the symmetry proper- 



ties of a system and its regularity or chaoticity. We show that the earlier suggested chaoticity 
criterion characterizes the initial symmetry breaking and destruction of the corresponding 
integrals of motion in a perturbed system, which leads to chaotisation. We compare our 
approach with known criterion of existence of the approximate quantum numbers by Hose 
and Taylor [5, 6] which is based on the analysis of effective Hamiltonians. 



$_i ■ One may ask if the new integrals of motion might appear in the perturbed system and 
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how can they influence the system's dynamics. The problem of construction of the ap- 
proximate integrals of motion has always attracted great attention due to its practical and 
philosophical importance. In classical mechanics we have KAM-theory which guarantees 
existence of invariant tori under perturbation; as to quantum mechanics the situation is 
somewhat tangled [5, 6]. The normal form method [7, 8] is well known for construction of 
the new approximate integrals of motion and integrable approximations to Hamiltonian. Its 
generalization on quantum systems was done in refs. [9-11]. After analyzing the reasons of 
divergence in the normal form method we propose a rather simple way for construction of 
approximate quantum numbers based on the unitary transformation of the basis integrals of 
motion. 

In the final part of the paper we discuss the possibility of the exact integrals existence 
in the perturbed system and formulate the hypothesis of formal integrability in stationary 
quantum mechanics. We define formal integrability as the existence (in mathematical sense) 
of complete set of independent mutually commuting operators which have some characteristic 
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differences from the usual first integrals of motion. It is impossible to write such operators 
in closed form (they have extremely complicated structure) and we can't find them before 
the solution of Schrodinger equation is obtained; these integrals don't correspond to known 
symmetry groups and are non-separable. We discuss why such a formal (mathematical) 
integrability doesn't contradict the system's quantum chaoticity. The conception of chaos 
means that the system has algorithmic complexity (it is very difficult to solve equations 
with high accuracy) and statistical hypothesis works well in the system. The existence of 
non-separable integrals of motion with very complicated structure can't help us to solve 
equations and doesn't influence the validity of statistical hypothesis. 

2 Destruction of quantum numbers and the chaoticity 
criterion 

Presently we consider stationary quantum system with Hamiltonian H as a sum of Hamil- 
tonian Hq of an integrable system and perturbation XV: 

H = H + XV. (1) 

Eigenfunctions {4> a } of the unperturbed integrable Hamiltonian Hq are common for some 
complete set of independent mutually commuting operators {J p }^ =1 (N - the number of 
degrees of freedom of the Hamiltonian H ). Eigenstates ipi of the full Hamiltonian may be 
expanded in eigenfunctions <p a of the unperturbed Hamiltonian Hq: 

^ = E»« = E c ^- ( 2 ) 

a a 

Let us consider the probability P a (Ei) to find the basis state <p a in the state ipi with energy 
which is equal to the squared absolute value of the corresponding coeflicient in (^) and 
define the energy width T" pr of P a (Ei) distribution: the minimal energy interval for which 
the sum of probabilities P a (Ei) is larger or equal to 0.5. Thus defined r" pr is the energy 
spreading width of basis state <j) a . The spectrum of Hq may be degenerate and then the 
irreducible representations of symmetry group of Hq consist of several basis functions which 
belong to one energy level (shell). We want to find a parameter characterizing the measure 
of initial symmetry breaking of Hq under the influence of the perturbation V. It's clear 
that such symmetry breaking results only due to significant mixing between functions from 
different irreducible representations. The mixing of states within one shell doesn't change 
their symmetry. Unless the spreading width T^ pr is smaller then the distance Dq between the 
neighboring levels in the spectrum of Hq we can distinguish "localization domain" (in energy) 
of one set of basis states from "localization domain" of another one. When the spreading 
width exceeds Dq we start loosing the " signatures" of basis functions in the spectrum of H 
and can't even approximately compare states ipi with irreducible representat ion of symmetry 
group Hq. Thus parameter 

ffi* = IVA) ( 3 ) 

is the natural measure of symmetry breaking. When the parameter a3 Q exceeds unity the 
symmetry of the Hamiltonian Hq disappears. Such a value of perturbation is accompanied by 
disappearance of the initial selection rules, the levels are distributed approximately uniformly 
(level repulsion) and the level spacing distribution approaches Wigner's law. One can say 
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that the transition from regularity to chaoticity has taken place in quantum system and se" 
may be considered as the parameter of chaoticity. 

The spreading width r" pr depends on the number a of basis state of Hamiltonian H Q , i.e. 
on its quantum numbers. The classical analogy to this is the dependence of the invariant 
torus stability on the corresponding vàlues of integrals of motion. It is clear that in order to 
obtain the global chaoticity characteristic in quantum case it is necessary to average over 
the basis states (pa belonging to the same irreducible representation (shell). The averaged 
chaoticity parameter r spr unlike local has an important feature of invariance with respect 
to the choice of the basis for integrable Hamiltonian Hq [12]. From the theoretical point of 
view the above chaoticity parameter has one more useful property. As it was shown in [3], 
in the semiclassical limit T^ w /% transforms into Lyapunov's exponent of the corresponding 
classical motion. 

Thus we see that criterion r spr "measures" destruction (fragmentation) of the irreducible 
representat ions of the basis and, hence of the Casimir operator (the main quantum number), 
which is the approximate integral of motion for small perturbations. The other quantum 
numbers of the basis states suffer destruction in general under a smaller perturbation due 
to the strong mixing within one shell. To define directly the degree of destruction of ap- 
proximate integrals of motion one may use an ordinary mean-square-root deviation from the 
mean value. This deviation of some operator A in the state ip a is calculated as follows: 



The operator A is the approximate integral of motion if the ratio of AA to the difference 
between its neighboring eigenvalues is less then 1. 

Another way to study and construct the approximate integrals of motion is the well known 
method of effective Hamiltonians. In the series of works Hose and Taylor [5, 6] suggested 
the criterion of existence of the effective Hamiltonian and of connected with it integrals of 
motion. According to this criterion we may build the convergent sequence of approximations 
to the effective Hamiltonian under the condition that projection of perturbed Hamiltonian 
wave function to the model space is greater then 0.5. Thus if projection of some states ip a of 
Hamiltonian H to the shell space exceeds 0.5, then the main quantum number has to be the 
approximate integral of motion for these states. It is obvious that this criterion practically 
coincides with our spreading width one. 

In order to compare the above criteria of destruction of the integrals of motion we have 
analysed the quantum Henon-Heiles system: 



The eigenfunctions were obtained using the oscillator basis of 496 states (30 shells). Fig.l 
shows the dependence on the perturbation intensity (energy E) of the parameter k, the 
averaged spreading width of operator N = n x + n y (see [|) and of the averaged projection 
(P s ) of the exact wave functions to the shell. For the sake of easier comparison with other 
quantities in Fig.l we plotted the value Pr = 2(1 — P s ). We see that destruction of the initial 
S , [/(2)-symmetry according to all these three criteria takes place approximately at E = 0.10. 

Thus T spr measures the degree of mixing between the irreducible representat ions and 
the destruction of the corresponding Cazimir operators. The question arises if some new 
integrals of motion might appear instead of the destroyed basis integrals. 




(4) 



H(q,p) = \{p\ + qï) + \{pl + qï) + Afe - ç 2 3 /3) 



(5) 
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The method of normal forms is the most known way to construct the approximate inte- 
grals and the integrable approximations to Hamiltonian in classical mechanics [7, 8]. It has 
been generalized on quantum systems in refs. [9-11]. In this method the perturbed wave 
functions are constructed as certain superposition of the basis functions belonging to a single 
irreducible representation. Therefore the symmetry of the wave functions is not changed and 
perturbation V with lower symmetry leads only to the splitting of the degenerated level. 

This approach gives rather good results when the perturbation is small. However Siegel 
(1941) proved the divergence of classical normal forms in the case of nonintegrable initial 
system. In quantum mechanics the question about convergence hasn't got a final solution, 
though the authors [10, 11] stress the asymptotic character of the series arising and confirm 
it by numerical calculations. 

Two reasons for the divergence in the method of normal forms may be pointed out. 
The first one is well known — it is the non-analyticity of solutions at the point À = 
(the replacement À — > — À in Hamiltonian may lead to significant changes of spectrum 
properties) and, as a consequence, to the divergence of the expansion into the powers of À. 
The second reason is as follows. In the quantum Birkgoff-Gustavson method the integrable 
approximations are constructed in the form of power series of operators which mix basis 
states only inside one irreducible representation of the basis symmetry group (one shell). If 
the interaction effectively mixes states from different irreducible representat ions the Birkgoff- 
Gustavson method obviously can't generate good integrable approximations and convergent 
integrals in principle. This is why the quantum numbers given by normal form might be 
good only when mixing between the different shells is weak (T spr < 1). 

The difnculties described doesn't mean the principal absence of approximate integrals of 
motion in the perturbed system, they only reveal the shortcomings of the methods. To get 
a convergent method of constructing the integrals one should use combinations of operators 
which mix states in any given finite-dimensional subspace. In this case we can improve the 
integrable approximation to H (in the sense of operator norm) by simply increasing the 
dimensions of this subspace. In the next section we'll describe the convergent procedure 
allowing to construct approximate integrals in a rather trivial way. 

3 The convergent method for integrals of motion 

Performing a unitary transformation U of some basis with a set of quantum numbers we 
can always easily find a new set of mutually commuting operators for which a transformed 
basis functions are eigenfunctions. One can see that J' = U .1 P U^ (J p - operators of initial 
basis) are the desirable operators and such bases are equivalent in the sense of the quantity 
of quantum numbers. If we assume the completeness of eigenstates ip a of Hamiltonian H 
in the Hilbert space H then ip a = U<p a (4> a is the eigenfunction of H ), because any two 
complete orthonormal bases are connected via some unitary transformation U . Operators 
J' p commute with H in the full Hilbert space TL. Actually, for any function from complete 
basis tp a 

[H, J' p ]4> a = HUJ p UH a - UJ p U^H^ a = HUJ p 4> a - E a UJ p <t> a = 

j pa HU(j) a - E a j pa U<f) a = j pa E a ip a - E a j pa ip a = 0. 

J' p are derived from J p with the aid of unitary transformation, hence they also form a 
complete set of independent commuting operators and their eigenvalues j pa uniquely define 
every eigenstate ip a of Hamiltonian H. Therefore H is a function depending on operators 
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J'\ H = H(J'), and E a = H(j pa ). Having determined the approximate wave functions 
of perturbed Hamiltonian we can construct the unitary transformation U and the approxi- 
mate integrals of motion J' . The question about convergence of the procedure suggested is 
reduced to investigating whether the corresponding method for solving of Schròdinger equa- 
tion converges or not. For example it has been proved that Ritz's method converges (with 
increasing basis) in the case of Hermitian operators with lower-bounded discrete spectrum 
[13]. 

The introduced operators J' seem to be formal unless we construct them explicitly as 
functions of dynamic variables (q,p). This, however is not difficult to do with the help of 
methods taken from the theory of continuous groups' representations. In the remaining part 
of this section the realization of this method is described in details and properties of the 
integrals obtained are discussed. 

Let us consider Schròdinger equation Hip = Eip with discrete spectrum E a and íp a and 
write the Hamiltonian H in the form of spectral decomposition 

ff = $> a VwC (6) 

a 

We represent the complete Hilbert space TC as a sum of the finite-dimensional model space 
V and its orthogonal adjunct Q: TC = V + Q. It is convenient to construct the space V 
of eigenfunctions {(j)^ d n=T by using some complete set of mutually commuting operators 
{ J p }p = i (N - the number of degrees of freedom H). We'll find the approximate wave functions 
of H in the P-space as a combination of basis states <p p (as, for example, in Ritz's variational 
method or in different versions of perturbation theory) 

VVa = E (7) 

Orthonormal states ip pa are derived from minimum condition for the energy functional in 
the P-space and they form a subspace in the P-space (obviously, only a small number of 
combinations (f7|) will satisfy Schròdinger equation with sufficient accuracy). We'll denote 
this subspace of solutions by S and the energy of states by E pa . The rest dimP — dimS of 
the basis functions in the P-space may be chosen arbitrarily, and we denote the new basis 
in the P-space {<^}^i P (0^ l = ip Pfí , ^ = 1, dimS). Let us show that the operator 

H s=Yl E P^ pa % a (8) 

(i) commutes with operators forming a complete set in the full Hilbert space 7í, i.e. it is 
integrable, (ii) H s is approximating H in the sense of operator norm in the iS-space, (iii) H s 
may be expressed in terms of dynamic variables (q,p) as well as the initial Hamiltonian H. 

By calculation of wave functions ippa we constructed simultaneously the unitary trans- 
formation cj)' = Ufifj, of space Ti, which is defined by coeficients in the P-space and is 
the identical transformation in the Q-space. Operators J' p = U J P U^ are known to form the 
complete set with the same quantum numbers j pa and eigenfunctions <f>' . As far as eigenfunc- 
tions of operators J' and H s in the 5-space coincide, the operators commute in this space. 
Outside the 5-space H s = and hence it also commute with J' . Therefore [H s , J'] = in 
the full space Tï. 

Now we are going to check that Hamiltonian H s is close to H in the sense of the operator 
norm in TC, i.. \\H — H s \\ s < e, under the condition that the residual of the approximate 
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solutions (H — E pa )ip pa = 5ip a doesn't exceed e: H^all < e. Really, for an arbitrary function 
X = Y l a 5 ïpp8, Hxll = 1, X^S 

\\Hx- H.x\ \ = 11(1] E a ip a ip* a - Epaippa^) <AM I = 

aeH aeS SeS 



1 1 a<5 ( E E a il) a ií)* a if) pS - E p sip p s) || <eJ2\ a5 \ < eV dimS. 
Ses ael·l Ses 

In the last estimate we used the fact that l a<5 | under the condition l a<5 | 2 — 1 reaches its 
maximum value when all as are identical and equal to 1/ \JdimS. The accuracy e depends 
on dimensionality of P-space; we may fix dimS and decrease e in such way that the norm 
||iï — H s \ \s should be as small as we need. 

The introduced operators H s , U and J' seem to be formal unless we construct them 
explicitly as functions of dynamic variables (q,p). Writing operator ippa^pa in terms of 
expansion (0) we have: 

H '= E {j2 E P^aC}^€- (9) 

n,ueV aeS 

The operator U may be also represented as a combination of basis operators (p^l'- 

U = l+ ^/^C+E^-i^r (io) 

a^pev iev 

The first dimS rows of the unitary matrix U a @ coincide with matrix c^, the rest dimV —dimS 
may be chosen in arbitrary way. It is easy to verify that components x don't change outside 
V when U acts on arbitrary state x = a^&s, while inside V they are transformed by unitary 
matrix: 

&0> aev f3ev 

Now we'll construct the basis operators as ordinary operators in the form of com- 
binations depending on variables (q, p) and acting on states of Hilbert space in co-ordinate 
representat ion. If Q is the group of transformations corresponding to the complete operator 
set J p {q,p), then P-space is in general a direct sum of irreducible representations T s of the 
group Ç: 

s 

Operator transforms the function <p u into (p^. If the group Q is Abelian one, the 

irreducible representations T s are one-dimensional and consist of the function <f> s . For non- 
Abelian group M and (fi u may belong to one irreducible representation. Our aim is to write 
the operators which generate all possible transitions between different T s and within some 
T s as well. It's not difficult to solve this problem by methods of group theory and actually 
the problem is equivalent to the realization of the basis M . 

To perform transformations inside the irreducible representations % it's suflicient to use 
combinations of generators of the basis symmetry group; to connect different % we need gen- 
erators of special non-invariance group of basis. Its infinite-dimensional irreducible represen- 
tation is spanned on our basis. We know non-invariance groups and corresponding algebras 
for various bases, for example, so(4, 2) — for Coulomb's basis and so(3, 2) — for isotropic 
two-dimensional oscillator basis [14]. After realization of operators of non-invariance àlgebra 
A a in the form of combinations of dynamic variables (q, p) we look for vacuum state 4>o for 
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which the decreasing operators from the set A a (q,p) give zero. The vacuum state forms 
one-dimensional irreducible representation of symmetry group Ç, and we shall naturally ob- 
tain states from other irreducible representations Q acting on it with creation operators from 
the set A a . Notation S^(A a ) defines the operator composed of generators A a which pro- 
duces basis function 0^: 0^ = S^fo S^cp^ = 0o- We don't present the general formula for 
= SUA a ) because it is not difficult to do it in any specific case (see section 4); usually 
are polynomials composed of generators A a the power of which increases with state's 
number /i. Then operator 0^0* on the Hilbert space 7í may be written as follows: 

^l = SlS u P v1 (11) 

where P v is the projector on the state <ft u . The projector P u may also be expressed in 
terms of dynamic variables (q,p) in the following way Let T(x) be the operators of unitary 
representation Q in Hilbert space, D s a Jx) — matrix elements of irreducible representation 
T s , dx — invariant Haar's measure on G. Then the projector P u on the basis state <p v G % 
may be presented as [15]: 

P v = dim% í dxD s *(x)T(x). (12) 

JG 

Operators (|T2|) are bounded, and since T(x) is the exponent to the power of generators Q 
which form sub-algebra with respect to operator àlgebra A a (q,p), we have achieved our goal 
- expressed the basis operators M 0*, H a (Q), U ([ÜJ) and integrals J' , in terms of variables 

Thus constructed integrable approximations H s (q,p) and integrals of motion J'(q,p), 
apart from their approximate character (the commutators with H are not exactly equal to 
zero) are local. The Hamiltonian H s is close to H in the sense of operator norm only in the 
finite-dimensional subspace S, while the operators J' are good invariants also only in the S- 
space ( outside the P-space they coincide with the old operators J p ). In the following section 
we shall demonstrate how the method works in Henon-Heiles problem and then continue to 
discuss the properties of the integrals obtained. 



4 Approximate integrals in Henon-Heiles problem 

Here we apply the method of the integral construction developed in the previous sections 
to the well known Henon-Heiles problem with Hamiltonian @. Introducing operators of 
creation and annihilation a\ = -^(qk + ipk)i a k — ^(^fc — Wk),k = 1,2 we construct, as 
usual, the Cartesian oscillator basis 

:(aln(4) n2 0o (13) 



We present the projector P ni „ 2 on the state ni r í2 (<?i5 Q2) = 0m(?i) ® 0n 2 (<?2) as a product of 
projectors on the states 4> ni (qi) and n2 (Ç2) of the corresponding one-dimensional oscillator. 
For one-dimensional oscillator Q = U(ï), T(x) = e m1ax , x e [0, 2ir], D s vv = e tux , where the 
number of state v equals the number of quanta in this state. Then according to (12]) 

Jo 2n 2w 

— (a)a — n)' 1 e™^ a ~ n ^ sin7r(a t a — n) = —(a)a — n) _1 sin7r(a t a — n). 
ix ix 



7 



We neglect the phase in the last expression because it does not affect the action of P n . It's 
easy to check that P n acts in the necessary way due to its slightly exòtic form: 



The total projector P nin2 takes the form 



Pnin 2 = -^{(A a i - n i) 1 (aÍ a 2-n 2 ) 1 sin 11 (a\a! - sm7í(a2a 2 - n 2 ). (14) 

7l z 



As a result the formulae (|ÏÏ|) together with (13) gives us the operator 



(a\ ) 111 ^ (a|)" 2 W (ai) ni M (a 2 ) ™ 2 M 
yjni(fj,y.n 2 (fj,)\ y / ni(z/)!n 2 (z/)! 



where ni and «2 are the quantum numbers of states. Using ([8]), (|Ï0| ) and determining 
coefficients and we get integrable approximation H s and approximate integrals J' p = 
U JpW . We may take n\ and n 2 or ni and n = ri\ + n 2 as independent integrals J p (p = 1,2). 
The coefficients of expansion and the energies E pa were calculated by Ritz's method with 
the aid of diagonalization of matrix H on the basis (|T3]). We shan't write explicit expressions 
for J' for they are very cumbersome: if D is the basis dimensionality, then operators J' 
consist of D 4 terms of the type (|Ï5D . To determine the degree of destruction of approximate 
integrals of motion J' in the states ip a we calculate the mean-square-root deviation J' with 
the help of @. The solutions obtained with the basis of 496 states (30 shells) were considered 
to be true wave functions. Fig.2 shows the dependence of averaged measure of destruction 
(fD of operators J'AJ\ = l J 2 = n) on the perturbation intensity (energy E) in the subspaces 
with different symmetries. Henon-Heiles Hamiltonian has a symmetry Cs v . Therefore the 
eigenfunctions' space can be divided into 4 subspaces A, B, C, D (the states belonging to C 
and D subspaces have the same energy and thus produce the sequences of degenerate levels). 
The approximate integrals J' were calculated for the P-space of different dimensionalities (1, 
10, 15 and 20 shells). One can see that increasing of P-space dimensionality is accompanied 
by the decrease of fragmentation of the approximate integrals in the iS-space; outside the 
iS-space operators J' loose their advantages in comparison with J p . 

This example shows that the approximate integrals J' really have smaller spreading 
(fragmentation) in the iS-space then the basis integrals of motion J p . We may get very small 
vàlues AJ' p w for the bounded states by increasing of dimV. As a result the analytical 
structure of J' becomes very complicated. 



5 Are integrability and chaos compatible? 

Now we consider the question about the convergence of the suggested procedure and the 
question about existence of exact integrals of motion of Hamiltonian H in the full Hilbert 
space (integrability of H). With the help of Ritz's method we can find in principle any 
finite number of states with any finite accuracy and thus construct integrable approximation 
H s to H and integrals J' in any finite subspace with any desirable accuracy. The question 
is: can we obtain the full infinite spectrum of H by tending V — » 7ï, because there is 
an effect of systematic "delay" of 5-space dimensionality with respect to dimensionality of 
the model P-space (only states far from the boundary of the approximate spectrum are 
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reasonably accurate in diagonalization) . In other words, does the sequence of Hamiltonians 
PHP converge to the initial Hamiltonian: PHP — > H while P — > 1. Physically it seems to 
be so but we are in a difíiculty to give rigorous mathematical proof, because the sequence 
PHP is not Caushy's sequence in the sense of operator norm in H. 

We can formulate the following hypothesis about integrability of H . If there exist a 
good (in mathematical sense) unitary operator U connecting two complete orthonormal 
bases: the initial cj) a and the basis of eigenstates ip a , then the Hamiltonian H, as we have 
seen above, commutes with the complete set of independent operators J' p and it may be 
expressed in terms of only these variables H = H(J' p ). (Moreover, according to Dirac [16] 
any functions of arbitrary complete orthonormal basis are eigenfunctons of some complete 
set of commuting observables. This allows to extend our conclusions for any observable, 
including the case of continuum spectra). Therefore the Hamiltonian H seems to be formally 
integrable. However, in the above considered example (the Henon-Heiles Hamiltonian) the 
system is chaotic according to all the criteria of Section 2. Moreover, we know for sure 
that in the classical limit this system is one of the textbook examples of chaoticity. The 
problem is how to remove the contradiction between the seeming formal integrability and 
the chaoticity of the system. 

One possible answer is connected with the properties of the new integrals J' p . These 
integrals of motion are independent and global (provided the convergence of PHP discussed 
above would be proved). However they have extremely complicated structure and can't be 
expressed in closed form. Therefore they are useless to separate variables (non-separable) 
and to solve a problem. We restore them after the approximate numerical solution has been 
found. These integrals don't give selection rules for transitions between levels. Therefore 
they are defmitely not the quantum analogs of the classical first (isolating) integrals, which 
define the classical regular integrable system. 

6 Conclusions 

The problem of existing and constructing of integrals of motion in stationary quantum 
mechanics and its connection with notion of quantum chaoticity has been investigated. It 
has been shown that the previously suggested quantum chaoticity criterion characterises 
destruction of initial symmetry of regular system and basis integrals of motion under the 
influence of perturbation. Our approach conforms with known probability criterion of Hose 
and Taylor [5-6] and direct estimate of fragmentation (f|). 

We use variational Ritz's method for explicit construction of approximate integrals of 
motion in the form of combinations depending on operators (q,p) though in principle an- 
other method for solving Schrodinger equation may be used. As a result we obtained finite 
large-dimensional sums consisting of non-invariance àlgebra operators in various powers and 
projectors nontrivially expressed in terms of invariance àlgebra generators. The quality of 
approximate integrals of motion is simply controlled by dimensionality of model space in 
use. 

These integrals of motion are independent and global (provided the convergence of 
PHP — > H discussed above would be proved). However they have extremely complicated 
structure but have extremely complicated structure and can't be expressed in closed form, 
therefore they are useless to separate variables (nonseparable). This also explains why the 
existence of these integrals doesn't create obstacles to statistical description of quantum 
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system. That's why such formal integrability (even if we'll prové the existence of global 
integrals rigorously) doesn't make system to be regular in the sense of absence of chaoti- 
cal properties. Therefore they are definitely not the quantum analogs of the classical flrst 
(isolating) integrals, which define the classical regular integrable system. 

One of the authors (IBI) is indebted to Prof.Zikiki and to the Organizing Committee of 
V.Gribov's Foundation for their support. 
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